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We demonstrate the implementation of Grover's quantum search algorithm on a liquid state 
nuclear magnetic resonance (NMR) quantum computer using essentially pure states. This was 
achieved using a two qubit device where the initial state is an essentially pure (e = 1.06±0.04) singlet 
nuclear spin state of a pair of H nuclei arising from a chemical reaction involving para-hydrogen. 
We have implemented Grover's search to find one of four inputs which satisfies a function. 

PACS numbers: 03.67.Lx, 82.56.-b, 03.67.Mn 



I. INTRODUCTION 

Quantum computers Q] are of enormous interest be- 
cause of their ability to perform calculations which ap- 
pear to be intractable on any conceivable classical com- 
puter. The leading technology for the implementation of 
quantum logic gates is liquid state nuclear magnetic res- 
onance (NMR) and NMR quantum computers 
HL IE 111 have been used to implement a range of quan- 
tum algorithms, mostly notably Shor's algorithm which 
was used to factor fifteen |9| . Unfortunately NMR devices 
have one great disadvantage: the small energy gap be- 
tween the Zeeman levels that are used as computational 
basis states means that direct cooling is not a practical 
method for obtaining a pure initial state. The most com- 
mon way round this is to prepare a pseudopure state [j| , 
whose behaviour is identical to that of a pure state up to 
a scaling factor, but this is not a good solution for two 
reasons. Firstly, the efficiency of pseudopure state prepa- 
ration falls off exponentially with the number of qubits 
in the quantum computer [lOj . which means that NMR 
devices cannot be scaled up to useful sizes, and secondly 
the states used are provably separable |n|, and so cannot 
exhibit the phenomenon of entanglement. 

These problems are not inherent in NMR quantum 
computation, but arise from the use of pseudopure states 
prepared from high temperature thermal states. We have 
recently shown how high purity initial states can be pre- 
pared using para- hydrogen techniques 0, 0] , and that 
these states can be used to perform the simplest quan- 
tum computation, Deutsch's algorithm |14| . Similar re- 
sults have been obtained by Hiibler et al. |l5j, albeit 
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with a lower spin state purity. Here we describe the use 
of a new molecular system which allows essentially pure 
initial states to be generated, and the implementation of 
Grover's quantum search algorithm on our pure state 
computer. State initialization on demand is achieved by 
using a 12 ns laser pulse to initiate a rapid chemical re- 
action, and is therefore both well controlled and time 
coherent with the NMR radiofrequency pulses. 

II. PAiM-HYDROGEN 

We prepare pure spin states using an effect call ed y ara- 
hydrogen induced polarization (PHIP) [3 El El HI • 
The existence of the para isomer of dihydrogen, H2, is 
a consequence of the Pauli principle, which requires the 
overall wavefunction to be antisymmetric with respect 
to exchange of the fermionic 1 H nuclei. It follows that 
H2 molecules in even rotational states, most notably the 
J = ground state, possess an antisymmetric nuclear 
spin wave function and correspond to nuclear spin sin- 
glets, termed para. Thus if H2 is cooled to a temperature 
of 20 K in the presence of a paramagnetic catalyst then 
essentially pure para-hydrogen will be obtained, with the 
singlet spin wavefunction 

hT>= ^(|oi>-|io», (1) 

where we have used the computational basis in which |0) 
corresponds to the ground state and |1) to the excited 
state of a spin-1/2 particle. The para-hydrogen molecule 
cannot be used directly for NMR quantum computing, 
due to its high symmetry, but this can be overcome by 
using a chemical reaction to prepare a new molecule, in 
which the two hydrogen atoms can be made distinct and 
can be separately addressed. For further details see [l2l 

13 El IHIII EI 13 Hi. 

In our previous work we have used the two hydride 
1 H nuclei in Ru(H) 2 (CO) 2 (dppe), where dppe indicates 
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FIG. 1: Quantum circuit to convert a singlet state jV^ - ) m t° 
the conventional initial state |00). H indicates a single qubit 
Hadamard gate, © indicates a single qubit NOT gate, and the 
initial two qubit gate is a controlled-NOT gate. 



l,2-bis(diphenylphosphino)ethane, as our NMR quantum 
computer. Detailed analysis of this system shows that 
the spins are not prepared in an absolutely pure singlet 
state, but in a slightly mixed state, which can be easily 
converted to a Werner state |2J of the form 



Pmt = (l~e)-+e\ifj ){i> |, 



(2) 



with a spin-state purity of e = 0.898 ± 0.026 PJ. Al- 
though this system is adequate for an initial proof of 
principle experiment, it is not ideally suited for quan- 
tum computing; in particular its H spin-spin relax- 
ation time (T2 = 0.58 s) is uncomfortably short. We 
have, therefore, prepared a novel but closely related 
system, Ru(H) 2 (CO) 2 (dpae), where dpae indicates 1,2- 
bis(diphenylarsino)ethane. This system has a slightly 
longer *H spin-spin relaxation time (T2 = 0.67s), and 
has the further advantage that the initial state can be 
prepared with a purity indistinguishable from one (the 
purity was measured using techniques described previ- 
ously mill and was found to be e = 1.06±0.04). Thus 
we describe our new system as a pure state NMR quan- 
tum computer. 

Although our quantum computer starts in a pure state, 
it does not start in the state conventionally assumed for 
two qubit quantum computers, that is |00) in the com- 
putational basis. It is, however, simple to convert our 
initial state to this form by using a disentangling circuit, 
such as that shown in Fig. ^ 



III. GROVER'S QUANTUM SEARCH 

Grover's quantum search algorithm [HI permits an ef- 
ficient search for one of k matching items within a search 
space of size n. The algorithm is most conveniently de- 
scribed in terms of a binary function, with n possible 
inputs of which k satisfy the function, that is produce an 
output of 1. The function is assumed to be implemented 
by means of an oracle, which will produce values of the 
function for any given input, but does not permit any 
other method of analyzing the function. The best clas- 
sical algorithm is simply to try inputs at random, and 
this will allow a satisfying input to be located in 0{n/k) 
queries. By contrast, Grover's quantum search allows a 
satisfying input to be located in 0(\Jn/k) queries. The 
simplest case occurs when k = 1 and n = 4, in which 



FIG. 2: A quantum circuit to implement Grover's search al- 
gorithm on a two qubit quantum computer. Following pre- 
vious practice [22| Hadamard gates H have been replaced by 
pseudo-Hadamard gates ft and ft— 1 , which correspond to 905- H 
rotations. The function being studied is encoded in the prop- 
agator Uf, and [Too replaces |00) by — 1 00) while leaving other 
basis states unchanged. At the end of the computation the 
two qubits are left in states corresponding to the satisfying 
input, that is f(pq) = 1. 



case the satisfying input can be located in a single query, 
and this is the case we will concentrate on here. 

The key element in an implementation of Grover's 
search is a quantum function evaluation oracle, Uf, which 
performs the transformation |22| 



\x) ^ (-!)/(»> \x) 



(3) 



where |x) indicates a quantum register in a state corre- 
sponding to the input x, and the function / has values 
of or 1. Typically x is written in binary form, and \x) 
is a set of qubits with values corresponding to successive 
bits of x. A quantum circuit which can be used to locate 
one satisfying input from four using this oracle is shown 
in Fig. H 



IV. THE EXPERIMENT 

Our two qubit system comprises the hydride H nu- 
clei in Ru(H) 2 (CO) 2 (dpae), where the hydride atoms are 
derived from para-hydrogen. The precursor compound 
Ru(CO) 3 (dpae) was prepared from Ru^ (CO) 12 and dpae 
using techniques described previously |12l |23(. Essen- 
tially pure para-hydrogen was prepared at a temperature 
of 18 K using a charcoal-based catalyst, and was intro- 
duced into a 5 mm NMR tube containing Ru(CO) 3 (dpae) 
dissolved in d6-benzene. The NMR tube was then trans- 
ferred into a 400 MHz spectrometer fitted with a 1 H/ 31 P 
probe modified for in situ photolysis |24| . The spectrom- 
eter triggered a 12 ns UV pulse of wavelength 308 nm 
from an MPB Technologies MSX-250 pulsed XeCl ex- 
cimer laser, irradiating the active region of the NMR sam- 
ple and producing the unstable species Ru(CO) 2 (dpae). 
Intermediates of this type have been shown to react 
with hydrogen on the sub-microsecond timescale pjj . 
which would in this case lead to the formation of 
Ru(H) 2 (CO) 2 (dpae). For further details see [H [HEi. 

The hydride resonances appear at — 7.61ppm and 
— 7.21ppm, with a frequency separation of 8 = 160 Hz, 
while the hydride J coupling ( 2 Jhh) was 4.8 Hz. The 
*H transmitter frequency was placed exactly between the 
two resonance frequencies. The laser flash acts as an 
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initialisation switch, generating the pure state on 
demand, which is subsequently used for the implementa- 
tion. 

The singlet state \ip~) can be converted to the desired 
initial state 1 00) using the NMR pulse sequence 
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P re P = [^] 90, [— ] 180, [— ] 180, [-] 90, (4) 

which has been described previously |14|. This pulse se- 
quence is composed of a series of hard RF pulses (which 
excite both spins equally) and periods of evolution un- 
der the background Hamiltonian, indicated by enclosing 
the appropriate evolution time in brackets. Following 
NMR conventions operations are applied from left to 
right. No selective pulses are used, but at two points 
equal and opposite z rotations are applied to the two 
qubits. These were achieved by evolution under the Zee- 
man Hamiltonian at frequencies of ±(5/2 Hz, with evo- 
lution under the small J coupling neglected over these 
short periods. NMR pulse sequences to implement the 
four possible functions / with a single satisfying value 
can be developed in the same way and are shown below, 
where each function is identified by listing its satisfying 
input. 



Poo 
Poi 
Pio 

Pll 



90_, 90, 90-, [i] 180, 



180, 
r 1 



180, 



i] [^] 180, [i] 180. 

— 1 180, [— ] 90_, 90, 90, 
4J j x l 4J j x y 



(5) 



Two of these functions require equal z rotations to be 
applied to the two qubits, and these were achieved using 
composite z pulses. 

The final element used in our pulse sequences is a gra- 
dient pulse, indicated by crush. This is a short period 
of evolution during which time the magnetic field is tem- 
porarily made spatially inhomogeneous, so that spins in 
different parts of the sample experience different Larmor 
frequencies. These crush pulses were applied immedi- 
ately before and after the Grover circuit, at which points 
the qubits should be in eigenstates of the computational 
basis. These states do not evolve under magnetic fields, 
and should not be affected by the crush pulse, but most 
error terms which can arise will be dephased [221 ]. 

Putting these elements together gives the final pulse 
sequence used to implement Grover's search, 



Pprep crush 90_, Pf 90, P( 



)0 90_, crush 90, acquire (6) 



where Pf corresponds to one of the four pulse sequences 
shown in equation [5] The state of the spin system is an- 
alyzed by applying a hard 90° pulse and observing the 
resulting NMR spectrum. With appropriate phasing the 
state of each qubit is then indicated by multiplets point - 



ref 




FIG. 3: Experimental spectra from our implementation of 
Grover's quantum search. The top spectrum is a phase ref- 
erence acquired from our computer in the state 1 00) . The 
four lower spectra were acquired after implementing Grover's 
search with the oracle set to each of the four functions in turn. 



ing upwards for |0) or downwards for |1). A phase refer- 
ence spectrum is easily obtained using the sequence 



Pprep crush 90, acquire 



(7) 



and the experimental results are shown in Fig. [21 

These spectra show exactly the pattern expected: each 
computation ends with the system in one of the four basis 
states, and the final state corresponds with the satisfy- 
ing input of the corresponding function. Note that in 
this paper we have used the same conventions as in |14| . 
where the first qubit occurs on the right hand side of the 
NMR spectrum; this is the reverse of the convention used 
in 22] . There are small imperfections visible in the spec- 
tra: in particular there are imbalances in the intensities 
of the two multiplets and of the two lines in each multi- 
plet, and the overall signal intensity is lower in spectra 
obtained from quantum computations than in the phase 
reference spectrum. The imbalances can be ascribed to 
errors in the implementations of the logic gates, while the 
signal loss is largely a consequence of relaxation during 
the pulse sequences. The effects of relaxation are less 
marked in this study than in our previous implementa- 
tion of Deutsch's algorithm [l4|; this partly reflects the 
longer relaxation times in our new spin system, but is 
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also a consequence of the symmetry in the gates applied 
to each qubit in Grover's algorithm. 

These results confirm that it is possible to use para- 
hydrogen techniques to implement quantum algorithms 
on a pure state NMR quantum computer. We are now 
seeking to extend these results to larger spin systems and 
more complex algorithms. 
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